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Abstract 

Understanding the biology and conducting effective conservation of migratory 
species requires an understanding of migratory connectivity - the geographic 
linkages of populations between stages of the annual cycle. Unfortunately, for 
most species, we are lacking such information. The North American Bird Ban- 
ding Laboratory (BBL) houses an extensive database of marking, recaptures and 
recoveries, and such data could provide migratory connectivity information for 
many species. To date, however, few species have been analyzed for migratory 
connectivity largely because heterogeneous re-encounter probabilities make 
interpretation problematic. We accounted for regional variation in re-encounter 
probabilities by borrowing information across species and by using effort cova- 
riates on recapture and recovery probabilities in a multistate capture-recapture 
and recovery model. The effort covariates were derived from recaptures and 
recoveries of species within the same regions. We estimated the migratory con- 
nectivity for three tern species breeding in North America and over-wintering 
in the tropics, common (Sterna hirundo), roseate (Sterna dougallii), and Cas- 
pian terns (Hydroprogne caspia). For western breeding terns, model-derived esti- 
mates of migratory connectivity differed considerably from those derived 
directly from the proportions of re-encounters. Conversely, for eastern breeding 
terns, estimates were merely refined by the inclusion of re-encounter probabili- 
ties. In general, eastern breeding terns were strongly connected to eastern South 
America, and western breeding terns were strongly linked to the more western 
parts of the nonbreeding range under both models. Through simulation, we 
found this approach is likely useful for many species in the BBL database, 
although precision improved with higher re-encounter probabilities and stron- 
ger migratory connectivity. We describe an approach to deal with the inherent 
biases in BBL banding and re-encounter data to demonstrate that this large 
dataset is a valuable source of information about the migratory connectivity of 
the birds of North America. 



Introduction 

Understanding the biology and conducting effective con- 
servation of migratory species requires knowledge of 
migratory connectivity, the geographic linkage of individ- 
uals, or populations between phases of the annual cycle 
(Marra et al. 2010). Quantifying the degree to which indi- 
viduals from a breeding population move to the same 
nonbreeding region is necessary for understanding how 
events during one phase of the annual cycle influence 
subsequent phases (Marra et al. 2005; Runge and Marra 



2005; Webster and Marra 2005). For example, habitat or 
precipitation experienced during the nonbreeding periods 
can influence an individual's survival to breeding (Wilson 
et al. 2011), arrival timing (Tottrup et al. 2008, 2012; 
McKellar et al. 2012), or reproductive success (Norris 
et al. 2004). Therefore, knowledge of migratory connec- 
tivity is essential for understanding population regulation, 
projecting the impacts of future threats, and solving diffi- 
cult environmental issues such as the spread of infectious 
disease (Webster and Marra 2005; Sheehy et al. 2010; 
Henkel et al. 2012). 
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Migratory Connectivity from Band Re-encounters 

Despite its importance, relatively little is known about 
the migratory connectivity of North American birds 
(Veen 2013). At present, we know the basic breeding and 
stationary nonbreeding ranges of most North American 
breeding species, but we do not know how populations 
are linked between these areas, including some of the 
most studied species in North America. Migratory con- 
nectivity has been quantified for a few North American 
species, particularly game species that migrate within 
North America (e.g., Diefenbach et al. 1988; Hestbeck 
et al. 1991). However, these studies typically provide esti- 
mates for only a few sites within the ranges of the species 
and have generally not accounted for the spatial variabil- 
ity in finding and reporting of re-encountered birds. 

New tracking devices are providing information about 
migratory connectivity (Hobson 2003; Wikelski et al. 
2007; Bridge et al. 2011; Douglas et al. 2012). However, 
current tracking methodologies have limitations. Satellite 
tracking devices remain too large for many species, and 
while they provide considerable data for each tagged indi- 
vidual, they are expensive to deploy on a large scale 
(Bridge et al. 2011). Geolocation dataloggers have limited 
precision and require that the animal be recaptured to 
retrieve data (Fudickar et al. 2011; Lisovski et al. 2012). 
Genetic markers, stable isotope signatures, morphology, 
and band recoveries are increasingly integrated to esti- 
mate migratory connectivity (Rundel et al. 2013; Rushing 
et al. 2014). These estimates from multiple data sources 
are more precise than estimates derived from single 
sources, but precision is still limited and modeling frame- 
works are not straightforward (Rundel et al. 2013). 

On the other hand, banding and re-encounter data are 
the most spatially accurate source of information on 
migratory connectivity for many species. The North 
American Bird Banding Laboratory (BBL) may be the 
largest inventory of tagged vertebrates in the western 
hemisphere (www.pwrc.usgs.gov/bbl). It houses over 
70 million banding and 4.5 million re-encounter records 
(live resightings or recaptures and dead recoveries) and is 
the only long-term dataset available for most Nearctic 
breeding bird species. The BBL database records include 
when, where, and how individuals were banded (from 
1955 to present) and re-encountered (from 1914 to pres- 
ent). The BBL database maintains the original banding 
records necessary to build individual capture histories, 
data not available within the European EURING databank 
(www.euring.org, Korner-Nievergelt et al. 2012). 

Given its availability as a source of information about 
migratory connectivity, banding data have too often been 
dismissed and underutilized (Korner-Nievergelt et al. 
2010a, 2012). Most efforts to interpret long-distance geo- 
graphic linkages from banding and re-encounter records 
have mapped the raw data (e.g., Brewer 2000; Bonlokke 
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et al. 2006; Sharrock 2010). However, these maps cannot 
be directly equated with population distributions because 
they do not account for the geographic variation in band- 
ing effort and re-encounter probabilities (Nichols 1996; 
Kendall et al. 2006; Korner-Nievergelt et al. 2010a). Spatial 
heterogeneity in finding and reporting of banded birds has 
long been documented (Nichols et al. 1995; Royle and 
Dubovsky 2001). For example, there is a high concentra- 
tion of common tern (Sterna hirundo) band recoveries in 
Guyana where they have been frequently trapped for food 
(Hays et al. 1997). However, common terns may be equally 
abundant in other parts of their nonbreeding range where 
trapping does not occur. Thus, interpretation of re- 
encounter locations requires accounting for the variation 
in the data due to the unknown observer and reporting 
distribution (Korner-Nievergelt et al. 2010a). 

There are well-developed statistical techniques for esti- 
mating movement probabilities from capture-recapture 
and recovery data that deal with the influence of the 
observer process by incorporating parameters for re- 
encounter probabilities (Brownie et al. 1993; Gimenez 
et al. 2007; Gauthier and Lebreton 2008). However, the 
small proportion of birds banded in breeding areas that 
are re-encountered in nonbreeding areas often limit the 
number of parameters that can be estimated with these 
data (Korner-Nievergelt et al. 2010a,b). Therefore, statisti- 
cal methods to overcome the issue of heterogeneous re- 
encounter probabilities in the face of small sample sizes 
are necessary to make large-scale banding and re-encoun- 
ter data an available source of information about migra- 
tory connectivity. 

In this paper, we present a novel and broadly applica- 
ble approach for deriving migratory connectivity estimates 
from banding and re-encounter data that accounts for 
heterogeneous re-encounter probabilities with limited 
data. We illustrate our approach with three species of tern 
that breed in North America, common tern, roseate tern 
(Sterna dougallii), and Caspian tern (Hydroprogne caspia). 
Finally, we use simulated data to evaluate the applicability 
of the method to species that vary in the number of indi- 
viduals banded and their re-encounter probabilities. 

Methods 

Model development 

The objective of this study was to estimate migratory con- 
nectivity from breeding to stationary nonbreeding regions 
using banding capture, recapture, and recovery data. We 
used Burnham's live-recapture dead-recovery modeling 
framework (Burnham 1993; Williams et al. 2002) which 
has four parameters, survival probability ((f)), recapture 
probability (p), recovery probability (r), and migration 
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(transition) probability (n, Table 1). The model incorpo- 
rated all individuals banded in breeding regions and any 
re-encounters in stationary nonbreeding regions. We ana- 
lyzed data in the R (R Development Core Team 2012) 
package RMark (Laake and Rexstad 2008) which is an 
interface for program MARK (White and Burnham 
1999). 

The state transition occurred once, from one breeding 
to one nonbreeding region. We model the encounter his- 
tories of birds captured and released and either recap- 
tured, resighted, recovered, or never re-encountered 
again. For example, the capture histories below are in the 
live dead (LD) format, consisting of three pairs of col- 
umns representing three capture occasions followed by 
their probability structure (Burnham 1993; White and 
Burnham 1999): 

L(A00002) = 7tA2<?> 1/2 (l - P 2 M1 -P2XI - q>)r 2 

L(C00030) = 7r C3 <p 1/2 (l -p 3 )q>p 3 [(l - <p)(l - r 3 ) + <p] 

The first occasion is for half of a year, from summer 
breeding to winter nonbreeding. Consequently, survival 
was calculated for half of the year for the first occasion, 
breeding to nonbreeding season, and was annual thereaf- 
ter, nonbreeding to nonbreeding season. During the first 
occasion, birds were banded in breeding areas A and C, 
respectively. The first dead occasion is always zero 
because recovery within the breeding season was not 
incorporated in the model. In both of the capture histo- 
ries, birds were not re-encountered live or dead during 
the first winter after banding. The bird in the first exam- 
ple banded in breeding region A was found dead in non- 
breeding region 2 during the second winter. The bird in 
the second example banded in breeding region C was 
recaptured or resighted in nonbreeding region 3 during 
the second winter. We followed the Burnham model for- 
mulation which makes it possible for a recapture and 
recovery to occur in the same year. For simplicity, the 
model assumes that recaptures happen during discrete 
capture occasions, but recoveries are not restricted to 
those occasions (Burnham 1993; Williams et al. 2002). 



This is generally not the case for large-scale banding data, 
where both recaptures and recoveries happen throughout 
the nonbreeding season. However, estimates of migratory 
connectivity were robust to this violation in preliminary 
simulations (J. Hostetler, unpublished analysis). There- 
fore, in the first capture history example, the bird survives 
one and one half years {(j) m 0(l-</>)) because recovery is 
assumed to have occurred between years. In the second 
example, the bird either dies after the last capture but is 
not recovered ((1— (f))(l— r 3 )) or survives ((/>). 

We increased parameter identifiability by estimating 
parameters associated with the sampling process, recap- 
ture (pj), and recovery (rj) probability, as the same among 
similar sized species that occupy overlapping nonbreeding 
habitats and regions (Thorup and Conn 2009). When one 
of these species did not occur in one or more of the non- 
breeding regions, we fixed migratory connectivity parame- 
ters to zero (White et al. 2006), providing known 
estimates for those 7t i; - (Thorup and Conn 2009; Korner- 
Nievergelt et al. 2010b). 

BBL re-encounter data for nongame species are not sys- 
tematically acquired; rather they are the long-term result 
of a combination of information from local-scale research 
projects, large-scale monitoring programs, and public 
reporting of band recoveries. Therefore, we estimated 
both recapture and recovery probabilities because the pro- 
cesses behind the recapture and recovery probabilities are 
likely to differ (Kendall et al. 2006). We took advantage 
of the large BBL database to build "effort" covariates for 
nonbreeding regions. The effort covariates borrowed data 
across many species to estimate the likelihood that a 
banded bird would be re-encountered in a nonbreeding 
region during the over-wintering period. We extracted 
records from the BBL database that occurred during the 
over-wintering months within the nonbreeding ranges of 
the species modeled. Of those records, we included spe- 
cies that were similar in size and habitat affiliation to the 
modeled species. The migratory connectivity of these spe- 
cies was not modeled due to low species-specific re- 
encounter numbers. The region-specific live (Pj) and dead 
(Ri) effort covariates (Table 1) were calculated as the total 



Table 1. Description of parameters used in the multistate live and dead re-encounter models. 
Name Description 

Kg Probability of a bird breeding in site / spends the nonbreeding season in site j 

Tj Probability that a bird dead in nonbreeding region j is found and reported 

Pj Probability that a bird alive in nonbreeding region ] is seen or captured and reported 

<f> Annual apparent survival probability 

Rj Site-specific effort covariate of dead re-encounter and reporting probability derived from 

the number of individuals of many species recovered in region j during nonbreeding 

Pj Site-specific effort covariate of live re-encounter and reporting probability derived from the number of individuals 

of many species recovered in region j during nonbreeding 
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number of re-encounters that fit our criteria in each 
region divided by the area of that region. They were 
incorporated into the re-encounter probabilities (pj and 
rj, respectively) with a logit link (Data SI). 

Application of this model requires assumptions impor- 
tant to consider when selecting species, and defining the 
spatial and temporal scales of interest. First, multistate 
mark-recapture models assume individuals within a state 
and capture occasion all have the same probability of 
detection, survival, and transition to the next occasion. 
The models also assume that the fate of each bird is cor- 
rectly assigned, with no tag loss, no influence from bands, 
and independence with the fate in the next occasion with 
respect to detection and survival. We also assumed that 
survival and detection during the nonbreeding period 
were independent of the breeding region. Re-encounter 
probabilities were assumed to be constant among species 
within over-wintering regions. Therefore, we chose species 
that overlapped in their nonbreeding ranges and were 
similar in size and habitat affiliations. We modeled spe- 
cies similar in size and habitat affiliations to minimize 
variability between species in finding and reporting of 
recoveries. 

The temporal scale selected was assumed to incorporate 
the full time periods when birds occurred in stationary 
regions, breeding and nonbreeding. We incorporated all 
available banding and re-encounter data, but reduced the 
number of parameters in our model by constraining them 
to be equal over years. Thus, we assumed that there was 
no temporal variability in migratory connectivity, sur- 
vival, or re-encounter probabilities from 1955 through 
2011. While re-encounter records are available as early as 
1914, computerized banding data are available and 
records were included for 1955 through 2011. Shorter 
time periods may be more appropriate for some BBL spe- 
cies (Visser et al. 2009). 

The spatial scales of nonbreeding regions were chosen 
to incorporate the full nonbreeding regions of species and 
individuals. We assumed that the entire nonbreeding 
range was incorporated to ensure that migratory connec- 
tivity parameters sum to one. We also assumed that indi- 
viduals in this model used only one breeding and one 
nonbreeding region during their lifetime. To avoid viola- 
tion of this assumption, we assigned large-scale breeding 
and nonbreeding regions, and evaluated these designa- 
tions with BBL data. No movement of any individuals of 
the study species was observed between breeding regions 
or between nonbreeding regions. 

Example using three tern species 

Common, roseate, and Caspian terns are long-distance 
migrants that breed in North America and have largely 
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overlapping nonbreeding ranges (Gochfeld et al. 1998; 
Cuthbert and Wires 1999; Nisbet 2002). We estimated 
migratory connectivity for these three species from breed- 
ing areas in the Northeast coast ("Eastern"), around the 
Great Lakes and along the St. Lawrence River ("Central"), 
and the interior West and Pacific coast ("Western") of 
North America to the coasts of four nonbreeding regions, 
(1) the Southern U.S. along the Gulf Coast and Florida 
and the Caribbean (GULF.CARIB), (2) eastern South 
America (ESAM), (3) Mexico and Central America 
(CAM), and (4) western South America (WSAM). 

Common terns have been extensively banded in their 
Atlantic Coast breeding colonies but less so in Central 
and Western parts of their North American breeding 
range (Nisbet 2002). Roseate terns are concentrated in a 
few major breeding colonies along the North Atlantic 
Coast of North America (Gochfeld et al. 1998). Less is 
known about the nonbreeding distribution of roseate 
terns, but they occur together with common terns in 
South America and have also been detected in the Carib- 
bean (Nisbet 1984; Gochfeld et al. 1998; Hays et al. 
1999). Caspian terns have a broad breeding distribution 
on coastlines, and inland lakes and rivers in North Amer- 
ica but are locally uncommon in most parts of their range 
(Cuthbert and Wires 1999; Wires and Cuthbert 2000; 
Morris et al. 2010). 

We incorporated all individuals that were banded dur- 
ing the breeding region and season, and re-encounter that 
occurred in the nonbreeding range and season. These spe- 
cies are primarily on breeding grounds during May 
through August. Less is known about arrival and depar- 
ture timing in nonbreeding areas, but they have been 
consistently observed from November through February 
in their nonbreeding ranges (Nisbet 1984; Cuthbert and 
Wires 1999; Hays et al. 1999; Nisbet et al. 2011). We 
included only re-encounters within 10 years of banding, 
encompassing the life span of most individuals for all 
three species (Cuthbert and Wires 1999; Hays et al. 
1999). We excluded re-encounters where the status of the 
bird or band was unknown after the re-encounter. Re- 
encounter records were categorized as (1) resighted or 
recaptured and released alive, (2) found dead or killed 
upon capture, or (3) removed from the marked popula- 
tion because the band was removed or the bird was held 
in captivity (Table 2). 

We built effort covariates of the number of individuals 
of similar species re-encountered within the nonbreeding 
regions during over-wintering months. Common, roseate, 
and Caspian terns are known to occur on beaches and 
other habitats with other species in their nonbreeding 
ranges (Blokpoel et al. 1984; Hays et al. 1999; Olmos 
2002; Bugoni and Vooren 2005). Thus, we included re- 
encounters of species represented in orders also commonly 
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Table 2. The number of terns banded during breeding and the percentage of those re-encountered during non-breeding. Of those re-encoun- 
tered, the status of the bird was live, dead, or unbanded after the re-encounter. Re-encounters of the tern species and the species in the effort 
covariate were obtained by non-scientists (shot, found dead, trapped) or scientists (recapture or resight). The effort covariate includes re-encoun- 
ters of coastal species in the non-breeding regions during the winter. 

Re-encounters 

Banded 





Total 


Total, % 


Status 






How obtained 




Dead, % 


Live, % 


Unbanded, % 


Non-scientific, % 


Scientific, % 


Common tern 


1,059,357 


0.09 


43.3 


42.5 


14.2 


76.1 


23.9 


Roseate tern 


104,204 


0.21 


19.6 


72.1 


8.2 


60.7 


39.3 


Caspian tern 


75,580 


0.16 


83.7 


9.8 


6.5 


94.3 


5.7 


Effort covariate 






49.7 


44.6 


5.7 


64.1 


35.9 



found in coastal habitats (Charadriiformes, Gaviiformes, 
Podicipediformes, Gmiformes, Pelecaniformes, Suliformes, 
Ciconiiformes, and Procellariiformes). We excluded An- 
seriiformes because it includes many game species and we 
expected the re-encounter probabilities from hunted spe- 
cies to be different from species that are not commonly 
hunted, including the three tern species. We compare re- 
encounter probability estimates and standard errors with 
and without the effort covariate. 

Simulation to assess bias and precision 

Simulated data were used to evaluate the broader applica- 
bility of this model. We simulated data that reflect available 
sample sizes for species in the North American BBL dataset. 
We have information about the number of re-encounters 
of 355 of the approximately 400 species of migratory North 
American breeding birds in the BBL database. Of those, 77 
species had >100 breeding to nonbreeding re-encounter 
records. The number of individuals re-encountered is a 
function of the number of birds banded during breeding 
and re-encounter probabilities during the nonbreeding 
period. The degree of migratory connectivity, the extent 
to which birds from the same breeding area migrate to 
the same nonbreeding area (Marra et al. 2010), may also 
influence the number of re-encounters. We assessed the 
magnitude and direction of bias when the strength of re- 
encounter probabilities or migratory connectivity varied. 
Where data were poor, very low re-encounter probability 
and weak connectivity, we expected greater bias. Therefore, 
the 27 scenarios varied in the number of birds banded (40, 
400, and 500 thousand), re-encounter probability (moder- 
ate, low, very low), and the strength of migratory connec- 
tivity (weak, moderate, and strong). 

We simulated migratory connectivity of one species 
from four breeding regions to four nonbreeding areas. 
We made the assumption that the number of individuals 
banded was the same in each breeding area. Re-encounter 
probabilities were the same in three of the four nonbreed- 
ing regions. Re-encounter probability was higher in the 



fourth nonbreeding region. The numbers of birds migrat- 
ing to each nonbreeding region varied with the strength 
of migratory connectivity to that region. We modeled 
both recapture and recovery with the same effort covari- 
ate. There were 100 replicates of 27 scenarios for which 
we calculated the mean number of individuals re-encoun- 
tered, coverage (proportion of estimates with confidence 
intervals that overlap the true value), root-mean-squared 
error (RMSE, a measure of the difference between actual 
and estimated values), and bias (difference between the 
mean of the estimates and the true value). 

Results 

We compared model results with and without the effort 
covariate. Recapture and recovery probability error was 
smaller in the model with the effort covariate, in all but 
one case (Fig. 1). Therefore, we present migratory connec- 
tivity estimates only from the model with the effort covari- 
ate but provide the code to run both models (Data SI). 
The ratio between the point estimates of recapture and 
recovery was small in all but three cases (range of values 
0.90-4.65). However, recapture and recovery probabilities 
in Central America were higher in the model that included 
effort (39.17 and 151.03 times higher, respectively) and 
recapture in the Gulf and Caribbean was moderately 
higher (13.26 times higher). The proportion of banded 
birds that were re-encountered varied among the species, 
but were fairly consistent within species (Table SI). 

When re-encounter probability was accounted for, wes- 
tern breeding terns had greater connectivity to the south- 
ern US (Table 3). Migratory connectivity estimates for 
eastern and central breeding terns were moderately differ- 
ent when re-encounter probability was accounted for 
(Table 3). Eastern breeding terns were strongly connected 
to eastern South America, and Western breeding birds 
were strongly tied to the western part of their nonbreed- 
ing range (Table 3, Fig. 2). In fact, we found no linkage 
between western breeding birds and eastern South Amer- 
ica. Central breeding common terns were more broadly 
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I 1 



PRSs 



Ss 



PRSs 



Ss 



Model 



Figure 1. Re-encounter probability estimates 
for common (Sterna hirundo), roseate (Sterna 
dougallii), and Caspian terns (Hydroprogne 
caspia) in four nonbreeding range regions (see 
Table 3). Estimates are from two models, with 
(PRSs) and without (Ss) an effort covariate on 
live, p, and dead, r, re-encounter probabilities 
(±95% CI). Parameters are assumed to be the 
same all three species in both models. 



connected to eastern and western nonbreeding regions 
(Fig. 2). The least precise estimates of migratory connec- 
tivity were for western breeding common terns where 
fewer individuals were banded. The percentage of birds 
banded within the breeding areas that were ever re- 
encountered was low for all three species (Table 2), but 
similar between breeding regions for each species (Table 
SI). Annual apparent survival was higher for roseate 
(0.79 ± 0.02) and similar for common (0.63 ± 0.01) and 
Caspian terns (0.67 ± 0.03). 

The simulations demonstrate that the accuracy and 
precision of migratory connectivity estimates generally 
increased with re-encounter probabilities and number of 
birds banded (Fig. 3, Table S2). Coverage was 85-100% 



across scenarios except when re-encounter probabilities 
and number of birds banded were lowest and when re- 
encounter probabilities and number of birds banded were 
highest (Table S2). The difference between actual and 
observed values (RMSE) was highest when re-encounter 
probabilities and number of birds banded were low 
(Table S2). However, when data were poor, with very low 
re-encounter probability and weak migratory connectivity, 
the model overestimated the effect of the effort covariate. 
Migratory connectivity estimates to nonbreeding areas 
with higher re-encounter probabilities were biased high 
while estimates to areas with lower re-encounter probabil- 
ities were biased low. There was little evidence of bias for 
all other cases (Table S2). 
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Table 3. Migratory connectivity estimates (SE) from the proportion of re-encounters in each nonbreeding region, and the multistate live and dead 
re-encounter model with the effort covariate. Breeding areas are in North America on the Northeast coast (Eastern), around the Great Lakes and 
along the St. Lawrence River (Central), and the interior West and Pacific coast (Western). Nonbreeding regions are along the coasts of the South- 
ern U.S. and the Caribbean (GULF.CARIB), eastern South America (ESAM), Mexico and Central America (CAM), and western South America 
(WSAM). 





Common tern 






Roseate tern 


Caspian tern 






Breeding 












Nonbreeding 


Eastern 


Central 


Western 


Eastern 


Central 


Western 


ESAM 
CAM 
WSAM 
GULF.CARIB 


Proportion of re-encounters 
0.95(0.01) 0.12(0.03) 
0.01 (0.003) 0.38 (0.0S) 
0.01 (0.004) 0.30 (0.0S) 
0.03 (0.006) 0.21 (0.04) 
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Figure 2. Migratory connectivity estimates for common (A, Sterna hirundo), roseate (B, Sterna dougallii), and Caspian terns (C, Hydroprogne 
caspia) breeding in North America. The width of the line is proportional to the strength of the connectivity. Estimates ± SE are shown. Less than 
1% of Eastern breeding common terns also migrated to Mexico and Central America and western South America. 



Discussion 

Most of the breeding birds of North America are repre- 
sented in what is likely the largest database of tagged ver- 
tebrates in the western hemisphere. We explore the use of 
these data to derive estimates of migratory connectivity 
from models that account for re-encounter probability. 
Using simulations, we show that the model is likely appli- 
cable for BBL species with >40,000 birds banded. 
Although, precision of estimates improve as the number 
of birds banded, re-encounter probability, and the 



strength of migratory connectivity increases. There was 
also little evidence of bias except for when migratory con- 
nectivity was weak and re-encounter probability was very 
low. Re-encounter probabilities are likely to vary with 
habitat affiliations and body size but, in general, they are 
likely to be relatively low in tropical nonbreeding regions. 
Low re-encounter probabilities will, therefore, likely limit 
the accuracy and precision of migratory connectivity esti- 
mates for many species in the BBL database. 

For western breeding terns, model-derived estimates of 
migratory connectivity differed considerably from those 
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Figure 3. Migratory connectivity estimates from simulated data. Birds from each of four breeding areas (A-D) migrate to each of four stationary 
nonbreeding areas (1^1). Mean estimates (maximum and minimum values) are from 100 replicates of 27 scenarios. Scenarios varied in the 
number of birds banded in each of four breeding areas (10, 100, and 500 thousand), the strength of migratory connectivity, and re-encounter 
probabilities. The strength of migratory connectivity is a 4 x 4 matrix with 16 values (weak all jr,j = 0.25, moderate 7r,j = 0.10, 0.15, 0.20, 0.55, 
and strong 7r,j = 0.05, 0.05, 0.15, 0.75). One migratory connectivity parameter from each breeding area (A-D) is calculated as one minus the sum 
of the other three. So, 12 of the 16 migratory connectivity parameters are estimated and presented here. In each scenario, re-encounter 
probability is higher in one nonbreeding area (1 is very low: 0.0015, low: 0.01, moderate: 0.08) and the same in the other nonbreeding areas (2- 
4 are very low: 0.0002, low: 0.002, moderate: 0.01). The solid lines indicate the true values. 



derived directly from the proportions of re-encounters. 
Conversely, for eastern breeding terns, estimates were 
merely refined by the inclusion of re-encounter probabili- 
ties. The proportion of banded birds that were re-encoun- 
tered (in any nonbreeding region) was fairly consistent 
between breeding regions within species. However, the 
number of birds banded in breeding regions varied con- 
siderably. Fewer birds were banded in the west and, there- 
fore, fewer were also re-encountered. If this concentration 
of banding effort is common for other species, then preci- 
sion of estimates are likely to reflect these geographic 
biases. However, we also found that most breeding 



regions were strongly connected to one or more non- 
breeding region. If this pattern of strong connectivity is 
common, it should increase our ability to estimate migra- 
tory connectivity from these data (Korner-Nievergelt et al. 
2012). 

Migratory connectivity information from additional 
data sources, such as tracking devices, will be useful for 
verifying results of these models. Five Eastern common 
terns tracked using geolocators spent the nonbreeding 
season along the eastern coast of South America with 
some individuals using more than one area during 
the same year but not moving between the regions we 
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designated (Nisbet et al. 2011). Further, re-encounter data 
prior to 1955 were not included in the model because the 
numbers of birds banded were unknown. However, it is 
possible to estimate the numbers of birds banded from 
breeding to breeding re-encounters (Korner-Nievergelt 
et al. 2012), making it possible to compare estimates from 
these data. We divided the nonbreeding ranges into a few 
large regions in an effort to meet the assumptions that all 
potential nonbreeding regions were incorporated and that 
individuals did not change breeding or nonbreeding 
regions during their lifetime. We did not detect any terns 
changing nonbreeding regions either within or between 
years. Fidelity to natal colony was high for Eastern breed- 
ing common and roseate terns (Spendelow et al. 1995; 
Nisbet and Spendelow 1999). Western Caspian terns 
moved between breeding colonies, but remained within 
the breeding region (Wires and Cuthbert 2000). 

The proportion of banded birds of each tern species that 
were re-encountered varied among the species despite their 
overlapping habitat associations and similarity in size. Re- 
encounter probabilities are likely to differ somewhat 
between species due to differences in such things as body 
size, habitat associations, behavior, and research effort. It is 
difficult to assess the extent to which variability among spe- 
cies is due to species-specific differences in re-encounter 
probabilities or the difference in the proportional linkages 
of species to nonbreeding regions that differ in re-encounter 
probabilities. Future work could compare estimates when 
re-encounter probabilities are modeled as proportional, 
instead of equal, between the three species (Korner- 
Nievergelt and Hofer 2009). The assumption that there was 
no difference in re-encounter probabilities between breed- 
ing regions was also difficult to evaluate because re- 
encounter probabilities differ between the nonbreeding 
regions and the species migrate to the nonbreeding regions 
in different proportions. We also could not assess whether 
re-encounter probabilities were equal within nonbreeding 
regions with our dataset. Korner-Nievergelt et al. (2012) 
found that violation of this assumption did not strongly 
bias parameter estimates for one species, but this result 
may not be generalizable to other species. 

Currently, we know year-round ranges for many species 
but not how migratory populations are connected 
throughout the annual cycle. This missing information is 
needed to understand the factors that regulate migratory 
populations and to forecast how climate change will affect 
the biology of migratory species. With over 70 million 
birds banded and over 1.2 million new records added 
every year, the methods presented here make the BBL 
database an important untapped resource for building a 
broader understanding of many of the migratory birds of 
North America. These data can provide baseline estimates 
of migratory connectivity that can be verified or 



improved with additional information from band sighting 
databases (e.g., www.mybandedtern.org, www.banded- 
birds.org, Kendall et al. 2006; Gratto-Trevor et al. 2012), 
stable isotopes, genetics, and tracking devices (Boulet 
et al. 2006; Ryder et al. 2011; Fraser et al. 2012; Macdon- 
ald et al. 2012). Large-scale banding and re-encounter 
data are available for many species, and it is possible to 
deal with the biases in these data to estimate migratory 
connectivity. 
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Table SI. The number of birds banded in each breeding 
region and re-encountered in each nonbreeding region, 
and the percent of birds banded that were encountered. 



Table S2. The model precision and bias from 100 repli- 
cates of simulated data. 

Data SI. R Code for MSLiveDead Mark models and Sim- 
ulation analyses in this paper. 
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